Ultrafast Insulator-Metal Phase Transition in VO2 Studied by Multiterahertz 

Spectroscopy 
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The ultrafast photoinduced insulator-metal transition in VO2 is studied at different temperatures 
and excitation fluences using multi-THz probe pulses. The spectrally resolved mid-infrared response 
allows us to trace separately the dynamics of lattice and electronic degrees of freedom with a time 
resolution of 40 fs. The critical fluence of the optical pump pulse which drives the system into a 
long-lived metallic state is found to increase with decreasing temperature. Under all measurement 
conditions we observe a modulation of the eigenfrequencies of the optical phonon modes induced 
by their anharmonic coupling to the coherent wave packet motion of V-V dimers at 6.1 THz. 
Furthermore, we find a weak quadratic coupling of the electronic response to the coherent dimer 
oscillation resulting in a modulation of the electronic conductivity at twice the frequency of the 
wave packet motion. The findings are discussed in the framework of a qualitative model based on 
an approximation of local photoexcitation of the vanadium dimers from the insulating state. 
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I. INTRODUCTION 



The insulator-metal transition in transition metal ox- 
ides, first reported by Morin,— is characterized by a sharp 
decrease in resistivity by several orders of magnitude 
when the compound is heated above a critical temper- 
ature Tc. Later a large number of compounds have been 
found to exhibit this type of phase transition. Exten- 
sive experimental information has been collected and the 
physical mechanisms responsible for the transition have 
been relatively well understood in many cases. ^ Neverthe- 
less, the origin of the phase transformation is still under 
debate in some systems. VO2 is a classical example that 
has attracted considerable interest for fundamental rea- 
sons as well as for possible applications, since its critical 
temperature — 340 K (67 °C) is located in the vicinity 
of ambient conditions. 

The change of the electronic properties in VO2 is ac- 
companied by a remarkable modification of the structure 
from the highly symmetric rutile to the low-temperature 
monoclinic phase where tilted V-V pairs are formedii The 
phase transition is of first order as confirmed by the di- 
vergence of the heat capacity and the presence of latent 
heat.— The high-temperature rutile phase of VO2 (space 
group D^^) can be visualized in terms of a body-centered 
tetragonal lattice formed by the vanadium atoms as de- 
picted in Fig. [lla), with each metal atom surrounded 
by an oxygen octahedron.— Each octahedron has virtu- 
ally cubic symmetry with very small orthorhombic dis- 
tortions. Symmetry lowering during the transition into 
the monoclinic phase (space group G\f^) leads to the dou- 
bling of the unit cell and the formation of V-V pairsS-iS as 
illustrated in Fig.[IJa). The soft mode of the rutile lattice 
corresponding to this structural transformation is located 
at the point R = (1/2,0,1/2) of the Brillouin 
Correspondingly, the monoclinic phase can be mapped 



onto the rutile lattice by two normal vibrational modes 
of vanadium atoms located at the F-point as a result of 
the Brillouin zone folding in the monoclinic phase^iS 

Working out the exact details of the electronic band 
structure of VO2 poses a true challenge to modern many- 
body theories. In particular, the correct prediction of the 
insulating gap remains a difficult task because the hierar- 
chy of the microscopic degrees of freedom (spin, charge, 
orbital and lattice) contributing to the phase transition is 
not clear. However, the phenomenological band scheme 
proposed by Goodenoug h^^d^ provides a quite realistic 
description by using simple arguments based on molecu- 
lar orbitals. It serves as a starting point for the discussion 
of the electronic properties of VO2 . 

In metallic VO2, hybridization between the oxygen 
2p and vanadium Sd orbitals leads to a- and 7r-type 
overlap. The filled cr and tt states will be primarily of 
O 2p character, whereas the corresponding antibonding 
bands are dominated by the V Sd orbitals [Fig.[T](b)] . The 
overlap of vanadium orbitals along the tetragonal c^' axis 
parallel to the vanadium chains results in the formation 
of a narrow band commonly labeled d\\. The Fermi level 
is located within all of the three t2g bands which are still 
partially overlapping despite being no longer degenerate 
as shown in Fig. [Ijb). 

Upon entering the low-temperature insulating phase, 
the experimentally observed structural changes bring the 
vanadium atoms closer to the apex of the barely modified 
oxygen cage, thus increasing the hybridization. Now only 
the d\\ band is occupied and the dimensionality of the 
electronic system is effectively reduced to unity. More- 
over, photoemission experiments report a splitting of the 
d\ I bandii'i^ by 2.5 eV that eventuates in the formation of 



0.7 eV between the lower d\\ band and 



a band gap of 0.6 ■ 

the Cg band and the transition to the insulating statCf^ 
Several theoretical models have been proposed based 
on Peierl a^^d'' and Mott-Hubbar d^^' scenarios. The 
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FIG. 1: (color online), (a) Vanadium sublattice of the rutile 
(left) and monoclinic (right) crystal structure of VO2. (b) 
Schematic band structure according to Refs. [l^ and [Tsl . The 
red arrow denotes the photoexcitation by an optical pump 
pulse with an energy of 1.5 eV. 

former supposes a lattice instability caused by electron- 
phonon interaction within a system of reduced dimen- 
sionality. The latter is based on strong electron-electron 
correlation effects. Recent ab initio calculations on 
the basis of cluster dynamical mean-field theory^S have 
demonstrated that the interplay of strong electronic in- 
teractions and the dimerization of vanadium atoms in the 
low-temperature phase play a crucial role in forming the 
energy gap. In this scenario which may be described as 
a " many-body Peierls insulator" , strong Coulomb corre- 
lations in VO2 lead to the formation of dynamical V-V 
singlet pairs that are necessary to trigger a Peierls tran- 
sition. The calculated electronic structure and optical 
conductivity in both phases of VO2 agree well with ex- 
perimental resultsi^ t^'^ 

The delicate balance of interactions drives transition 
metal oxides, such as VO2, into a critical regime that is 
ruled by phase competition and that reacts exceedingly 
sensitively to external stimuli. When ultrafast photoex- 
citation favors one of the competing phases via the in- 
teraction of a photoexcited state with lattice, spin, or 
charge degrees of freedom, a dramatic phase conversion 
may occur. These phenomena are highly cooperative so 
that the structural relaxation processes of the electronic 
excited states are not independent, as in conventional di- 
lute excitonic or photo-chemical absorption, but entail a 



photoinduced phase transformation toward a new lattice 
structure and electronic order. ^"^-^^ This scenario opens 
the way for a light pulse to induce symmetry breaking 
from a stable phase and so to establish a new long-range 
order. In this context, ultrafast time-resolved techniques 
are emerging as a tool to study phase transitions in com- 
plex materialsi^ Time is introduced as an additional 
parameter. It promotes the possibility to unravel the 
contributing degrees of freedom and advance our under- 
standing of sophisticated physical processes. Evidently, 
the time-resolved technique of choice should be sensitive 
to the relevant degrees of freedom. 

The interest in the insulator-metal transition of VO2 
was rejuvenated after Becker et ali^ discovered that the 
metallic state of VO2 may not only be induced thermally 
but also optically on an ultrafast timescale. Later, X- 
ray diffraction^ and optical^ pump-probe experiments 
proved the non-thermal character of the transition and 
suggested the limiting time for switching to be less than 
100 fs. A large transient pump-induced increase of the 
conductivity in the THz frequency range was recently re- 
ported by several groups<^— Therefore, nowadays VO2 
has attracted considerable attention due to potential 
applications such as ultrafast control of light in VO2- 
based photonic crystals, ■^^ optical switches^"— and op- 
tical storage devices. '^^ Efficient harvesting of this large 
technological potential demands a thorough understand- 
ing of the microscopic physical processes. 

Recently, we reported measurements of the insulator- 
metal transition in VO2 induced by 12-fs optical pulses 
and probed by ultrabroadband multi-THz transients)^ 
These experiments made it possible to discriminate spec- 
trally between the excitation dynamics of electronic and 
lattice degrees of freedom, revealing their fundamentally 
different character. Based on these observations a novel 
qualitative model of the photoinduced insulator-metal 
transition in VO2 has been suggested. 

The present paper extends the previous short report 
by addressing such important aspects as a temperature 
dependence of the photoinduced changes and a detailed 
insight into the coherent oscillation dynamics. In this pa- 
per, we present a systematic study of the photoinduced 
transition at different temperatures and excitation flu- 
ences. The minimum energy of the pump pulse required 
to induce the cooperative non-thermal transition into the 
metallic state is found to be comparable to the thermo- 
dynamic energy difference between the monoclinic and 
rutile phases. This fact indicates an ultrafast switching 
of the lattice structure that instantaneously follows the 
electronic excitation of vanadium dimers. As a charac- 
teristic fingerprint of this process a modulation of the 
eigenfrequencies of the optical phonon modes induced by 
their anharmonic coupling to the coherent wave packet 
motion of V-V dimers is observed under all measurement 
conditions. In addition, a modulation of the electronic 
response at twice the frequency of the wave packet mo- 
tion gives evidence of a quadratic coupling to the struc- 
tural order parameter. The observed phenomena allow a 
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FIG. 2: (color online). Schematic of the optical pump / multi- 
THz probe setup. The output of the amplified Ti:sapphire 
laser system is split into three branches: (1) The optical pump 
pulse with variable delay r; (2) the THz probe pulse with 
variable delay T; and (3) the gating pulse for field-resolved 
electro-optic detection. The sample is mounted in a liquid 
helium cryostat equipped with CVD diamond windows. 



deepened insight into the interplay of electronic and lat- 
tice degrees of freedom during the insulator-metal phase 
transition in VO2. The qualitative model of the ultrafast 
phase transition proposed previously is discussed in view 
of the additional results. 



sample which is held at a preset substrate temperature 
Tl. The transmitted radiation is refocused by another 
set of two parabolic mirrors onto the 50 ^m GaSe electro- 
optic sensor. By varying the delay T of the gating pulse 
(branch 3) with respect to the transient, the time evolu- 
tion of electric field can be traced for selected delay times 
T of the pump pulse. Fourier analysis of the time-domain 
data yields two-dimensional (2D) maps of both absolute 
amplitude and phase of the pump-induced transmission 
change. The time resolution along the r-axis in our setup 
is determined by the detector bandwidth and set to be 
approximately 40 fs. The real and imaginary parts of the 
refractive index of the VO2 thin film are extracted from 
the measured equilibrium transmission and the pump- 
induced transmission change by numerical solution of a 
corresponding Fresnel equation,^! Any other optical con- 
stants, e.g. optical conductivity cr(w, r) or permittivity 
e(a;,T), are available from these data through basic elec- 
trodynamic relations 42. 

The spectrally-integrated dynamics of the optical con- 
ductivity is measured if the delay T is fixed at the maxi- 
mum of the multi-THz transient and only the delay of the 
pump pulse T is scanned. Since the total spectral inten- 
sity of the probe pulse is proportional to the peak value 
of the electric field, a measured pump-induced change of 
its amplitude characterizes changes in absorption or op- 
tical conductivity of a sample in the whole spectral range 
covered by the multi-THz transient. 



II. EXPERIMENT 



III. RESULTS 



The sample is a 120 nm thin film of polycrystalline 
VO2 grown by pulsed laser deposition on a CVD dia- 
mond substrate^i A vanadium-metal target is ablated 
in 250 mTorr of oxygen by a KrF excimer laser (repeti- 
tion rate 25 Hz) at a wavelength of 248 nm and a fluence 
of 4 mJ /cm^ . The as-grown films are non-stoichiometric 
with an approximate composition close to VO1.7. There- 
fore, the film is subsequently oxidized at 450 °C under 
7x 10"'^ mbar of oxygen pressure. After oxidizing to VO2, 
the stoichiometry is verified via Rutherford backscatter- 
ing and observation of electrical and optical switching, 
while crystallinity and phase homogeneity are confirmed 
by x-ray diffraction. 

Our multi-THz setup is based on a home-built low- 
noise Ti:sapphire amplifier system for intense 12-fs light 
pulses centered at a photon energy of 1.55 eV (800 nm 
wavelength). ■^^ Pulse energies up to 0.5 /iJ are generated 
at variable repetition rates of up to 4 MHz. The sim- 
plified schematic of the optical pump / multi-THz probe 
setup is depicted in Fig. O Part of the laser output de- 
noted as branch 1 is used for optical pumping of the 
sample. The beam in branch 2 generates phase-locked 
multi-THz probe transients by optical rectification in a 
50 /im thick GaSe crystal. '^^ A pair of parabolic mirrors 
focus the THz beam to the spot excited by the pump 
beam. The multi-THz transients are focused onto the 



A. Equilibrium optical response of the VO2 sample 

When the sample temperature is varied across the 
phase transition, the transmission in the THz range 
shows a pronounced hysteresis illustrated in Fig. [2Ia). 
The insulator-metal transition occurs at Tl = 339. 5±1 K 
upon heating and at 333 ± 1 K when cooling down. The 
resulting narrow hysteresis width of 6.5 K attests to the 
high quality of the thin film sample^^ The hysteretic 
behavior confirms the first-order character of the phase 
transition in VO2. Previous calorimetric studies demon- 
strated the presence of a latent heat of 241 3/cm^^ at the 
transition point All measurements presented in this pa- 
per were acquired on the heating part of the hysteresis 
cycle. 

Fig. ^h) shows the real part of the optical conduc- 
tivity cri(a;) for two selected sample temperatures be- 
low and above Tc. Although ai of insulating VO2 is 
expected to vanish, at Tl — 295 K it exhibits pro- 
nounced maxima at hw — 50, 62, and 74 meV (13, 15, 
and 18 THz). These peaks correspond to the transverse 
optical phonon resonances of the monoclinic lattice re- 
lated to vibrations of the oxygen cages surrounding the V 
atomSfS-d^ Vanadium-dominated normal modes, in con- 
trast, are known to oscillate in the low-frequency regime 
between 2 and 6 THz.^'^^ The spectral region above 
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FIG. 3: (color online), (a) Temperature dependence of the 
normalized transmittance through a VO2 thin film upon heat- 
ing (red guide to the eye) and cooling (blue guide to the eye) . 
(b) Real part of the optical conductivity of the VO2 sample 
below and above Tc. 



85 meV is free of infrared-active resonances and displays 
a low conductivity in the insulating phase. The spectral 
positions and weights of the phonon resonances in the 
monoclinic phase are in good agreement with previous 
infrared studies on VO2 single crystals. ^'"'^ The observed 
phonon resonances can be identified by comparing their 
eigenfrequencies and oscillator strengths with data re- 
ported by Barker et al4^ Two strong resonances at 62 
and 74 meV are assigned to the oxygen vibrations per- 
pendicular to the monoclinic um axis (collinear with cr of 
the rutile phase). The weaker resonance around 50 meV 
is related to the vibration polarized along the am axis. 
As expected for polycrystalline thin films, the phonon 
linewidths of about 10 meV observed in our sample are 
about twice larger than those of single crystalline samples 
(4-5 meV).— First, our thin-film sample are polycrys- 
talline and have reduced long-range order which causes 
a broadening of the observed lattice modes. Second, the 
experiments performed with polycrystalline thin films al- 
ways average over different crystal orientations. 

The optical conductivity in the metallic phase in- 
creases monotonically with temperature. However, it 
does not increase with decreasing energy as expected for 
a Drude response of conventional metals [Fig.jSfb)]. The 
peculiar spectral shape arises from the co-existence of 
insulating and metallic domains typical of a first-order 
transition. Hence, in a DC resistivity measurement, a 
sample appears metallic as soon as a percolation path 
connecting metallic domains is formed. On the other 
hand, the optical response averaged over many domains 
is determined by an effective dielectric function that de- 
pends on the volume fractions of the monoclinic and ru- 
tile phase. A good description of experimental infrared 
spectra using Bruggeman's effective medium theory^ has 
been demonstrated!^""— Based on this success, recent 
studies rely on the effective medium theory for a quanti- 
tative analysis of the optical spectrai ^^'^^" — 
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FIG. 4: (color online), (a) Spectrally integrated transient 
change of the THz conductivity after excitation by a 12-fs 
near-infrared laser pulse at Tl = 295 K for various pump 
fluences. (b) Fluence dependence of A(Ti(r) at r = 60 fs 
(blue crosses) and 1 ps (red triangles), (c) Extrapolation 
of A(Ti(l ps) curves (red triangles: 295 K, magenta circles: 
320 K) to a critical fluence of $^(295 K) = 4.6 mj/cm^ and 
"I'c(320 K) = 3.5 mJ/cm^, respectively, (d) Dependence of 
threshold fluence $c on lattice temperature Tl- The solid line 
is the thermodynamic energy difference between the metallic 
state and the insulating state at given temperature calcu- 
lated according to Ref. 3. The dashed line marks the phase- 
transition temperature Tc. 



B. Time-resolved spectrally integrated multi-THz 
response of photoexcited VO2 

Now we consider the non-equilibrium response of VO2 
induced by the optical pump pulse. The pump photons 
at 1.55 eV promote electrons from the spUt-off d\\ band 
below the Fermi level to the conduction band tt states 
[see Fig. [I{b)]. The optically generated free electron- 
hole pairs, in turn, give rise to a finite photoconductivity. 
Fig. Ufa) depicts the ultrafast dynamics of the spectrally 
integrated conductivity change A(Ti(t), for a series of 
pump fluences $ recorded at a sample temperature of 
Tl = 295 K. The initial increase of the infrared conduc- 
tivity peaks within a resolution-limited time interval of 
60 fs, marking the completion of the excitation process. 
At low excitation densities, a sharp onset of the pump- 
induced signal is followed by a non-exponential sub-ps 
decay. At higher fluences an increasing background of 
long-lived conductivity appears that is constant within 
a time window of 10 ps. The density of directly excited 
photocarriers contributing to the initial conductivity de- 
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FIG. 5: (color online). Decomposition of the decay dynamics, 
(a) Spectrally integrated dynamics of the THz conductivity 
after excitation at fluence "1? = 14 mj/cm^ and temperature 
Tl = 4 K (black circles). Subtracting the biexponential de- 
cay component (red line) reveals a coherent oscillation (blue 
circles) to which a cosine function with a frequency of [j/27r 
= 6.1 THz (green line) was fitted, (b) Spectrum of the oscil- 
latory component (blue circles) and fit by a Lorentzian func- 
tion (blue line), (c) The spectrum of the oscillatory modu- 
lation observed in a degenerate pump-probe experiment^^ at 
1.55 eV (red line) and the unpolarized spontaneous Raman 
spectrum^ (gray line). 



pends linearly on the fluence. Thus, we expect the quasi- 
instantaneous signal to scale linearly with $, as well. 
Extracting the values of the pump-induced conductiv- 
ity at T = 60 fs confirms this expectation, as shown in 
Fig. m^b). In contrast, the corresponding value at r = 
1 ps vanishes for small fluences while it increases rapidly 
above a threshold $c(295 K) = 4.6 mJ/cm^, which is 
determined by extrapolating the linear part of the latter 
graph to zero [inset of Fig.|4l[b)]. The threshold in excita- 
tion density separates two regimes: below $c tbe directly 
photo-generated electron-hole pairs populate delocalized 
states and thus give rise to the initial photoconductivity. 
However, the lifetime of these states is limited by a very 
effective trapping or relaxation mechanism that results 
in the rapid sub-ps decay of the quasi-instantaneous sig- 
nal. Increasing the excitation density above the thresh- 
old of only one absorbed photon per approximately 10 
V-V dimers triggers a cooperative transition from the 
insulating to the metallic state that renders the relax- 
ation pathway inoperative. The subsequent persistence 
of the photoconductivity indicates the transition to the 
metallic phase. However, the initial insulating state is 
restored before the arrival of the next pump pulse due to 
heat dissipation to the diamond substrate (the substrate 
temperature Tl was kept below the hysteresis region). 

Raising the sample temperature closer to Tc leads to 
a decrease of the threshold fluence $c [see Fig. Sfc)] and 
the photoinduced conductivity becomes more long-lived 
as compared to room temperature. On the other hand, 
the cooperative transition may be suppressed when the 



lattice is cooled to cryogenic temperatures although the 
excitation density exceeds its critical room-temperature 
value [see Sec. IIIIC] . Indeed, as it is shown in Fig. Wid) 
the fluence threshold experiences a significant reduction 
as the critical temperature is approached from below. 
Since the threshold fluence corresponds to the minimal 
energy of the pump pulse necessary to induce the coop- 
erative transition into the metallic phase, it is instruc- 
tive to compare it to the difference of the thermal energy 
of the VO2 lattice between the metallic and insulating 
phases. This difference can be calculated as the integral 
of the lattice heat capacity between a given temperature 
and Tc plus the latent heat of the insulator-metal transi- 
tion. The resulting curve calculated using the data from 
Ref. '3' is shown in Fig. lU^d) as a solid line. In order to en- 
able a direct comparison with pump fluence, the thermal 
energy density is recalculated in terms of a surface en- 
ergy density assuming a 100-nm-thick excited layer. The 
agreement between both data sets is remarkable. 

Since the electronic contribution to the speciflc heat 
in the studied temperature range is negligible compared 
to the lattice contribution, we can conclude that optical 
pumping above ^dT) should induce a structural trans- 
formation of VO2 into the rutile phase on the time scale 
shorter than 1 ps after the excitation which we have cho- 
sen for the estimate of the excitation threshold. Indeed, 
the dynamics of the optical conductivity in Fig.^l^a) con- 
firms this assumption. 

It should be noted that the agreement between optical 
pump and thermal energies [Fig. Iljc)] does not imply a 
thermal character of the photoinduced transition in VO2. 
Clearly, the sub-ps switching time would be untypically 
short for a thermally driven process. As a matter of fact, 
both coherent (optical switching) and incoherent (ther- 
mal heating) mechanisms are observed, but they are well 
separated in time. For example, after photoexcitation at 
295 K the pump-probe signal starts rising again after sev- 
eral tens of ps and reaches its maximum level after 300 ps, 
independent of the excitation fluence [not shown]. The 
phase growth in the thin fllm is driven by hot phonons 
propagating at the speed of sound and proceeds inco- 
herently at spatially separated sites. The conductivity 
dynamics of this thermal phase transition activated by 
laser heating was studied in detail by Hilton et al.^^ and 
is not a subject of this paper. 

A closer look at the decay dynamics of the THz con- 
ductivity reveals a fast oscillation superimposed on the 
monotonically decreasing signal [Fig.lSJa)]. The presence 
of this feature in the pump-probe signal is a flrst sub- 
stantial signature of the participation of lattice modes in 
the phase transition. Fig. [Sfa) shows how subtraction 
of a biexponential decay from the pump-probe signal re- 
veals a coherent cosine-like oscillation with a frequency 
of 6.1 THz. Apparently, the flrst cycle of the observed 
coherent oscillation is notably stronger and cannot be 
well described by a simple cosine function. This points 
out a strongly non-equilibrium character of the multi- 
THz response during the flrst oscillation period of 160 fs 
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FIG. 6: (color online). 2D optical pump / multi-THz probe 
data: (a) Equilibrium conductivity of insulating VO2 at 
295 K. Color plots of the pump-induced changes of the con- 
ductivity A(Ji(a;,T): (b) at Tl = 4 K and an incident finance 
of $ = 7.5 mJ/cm^ at Tl = 250 K and pump fluence (c) <E> — 
3 mJ/cm^ and (d) $ — 7.5 mj/cm^. The broken vertical lines 
indicate the frequency positions of cross sections reproduced 
in Fig. [71 



when the lattice distortion is accompanied by a strong 
change in electronic structure. The following harmonic 
oscillations describe a pure lattice vibration around a new 
quasi-equilibrium configuration. It should be noted that 
coherent oscillations at a constant frequency of 6.1 THz 
are observed at all lattice temperatures below Tc inde- 
pendent of the excitation fluence. Pumping the metallic 
phase {T > Tc) reduces the conductivity and the sign of 
the spectrally integrated response is inverted. Thus, for 
delay times r > 100 fs, the observed behavior is typical 
of a hot electron gas. 

The Fourier transform of the oscillatory component ex- 
cluding the first cycle is shown in Fig.[5l[b). It is well fit- 
ted by a single Lorentzian term. In contrast, all-optical 
ultrafast reflectivity experiments^^ have revealed impul- 
sive excitation of two vibrational modes at 5.85 THz 
and 6.75 THz which coincide with totally symmetric Ag 
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FIG. 7: (color online). Cross sections of Fig. [6] along the time 
axis r for a photon energy of (a) Tiuj — 60 meV and (b) hcj = 
100 meV. The curves taken at $ = 3 mj/cm^ are scaled up by 
a factor of 2.4. (c) and (d) The oscillating components of the 
cross sections through the 2D scan in Fig. [SJd) for a photon 
energy of (c) huj — 60 meV and (d) huu = 92 meV. Green solid 
lines: (c) fit of the oscillating component by a cosine function 
with frequency of 6.1 THz, (d) the fitting function shown in 
panel (c) squared with subtracted constant background. The 
oscillation occurs at doubled frequency of 12.2 THz. 



Raman modes of the monoclinic lattic o^^'^^ [Fig. [5l[c)]. 
These phonon modes are critical to the metal-insulator 
transition in VO2 and describe stretching and tilting of 
vanadium dimers which map the monoclinic onto the ru- 
tile lattice. Thus, the corresponding resonances are only 
observed in the ground insulating state of VO2 probed 
by spontaneous Raman spectroscopy. The observation 
of a single oscillation frequency instead of the two fre- 
quencies in the dimerized monoclinic phase indicates the 
higher lattice symmetry of the excited state probed in 
our experiment. 

In general, pump-probe techniques utilizing an impul- 
sive excitation may be sensitive to ground-state as well as 
to excited-state vibrationsi^^i^ This can be clearly seen 
in Fig. ^c) where the experimental data from a degen- 
erate femtosecond pump-probe experiment at a photon 
energy of 1.55 eV are shown. Besides two Ag Raman 
modes of the ground state of VO2 an additional shoulder 
appears between both maxima, exactly at the frequency 
of 6.1 THz where the coherent oscillation in the multi- 
THz probe occurs (see vertical dashed line in Fig. [SJ. 
In contrast, the multi-THz probe obviously is not sensi- 
tive to the coherent wave packet dynamics in the ground 
state of VO2 and probes almost exclusively the coherent 
oscillation in the excited electronic state. The reason for 
such selectivity cannot be understood based solely on the 
spectrally integrated data which do not unravel the mul- 
tiple participating degrees of freedom. This limitation is 
overcome in a full 2D optical pump / multi-THz probe 
experiment. 



7 



C. Spectrally resolved dynamics 

The spectrally resolved pump-probe response is ob- 
tained by measuring the complete time profiles of multi- 
THz transients as well as their pump-induced changes 
at varying pump-probe delay times t and subsequent 
Fourier analysis as described in Sec. [Ill The resulting 
2D plots shown in Fig. [S] depict the spectral profiles of 
the conductivity changes A(Ti(a;, r) as a function of time 
delay r at different temperatures and pump fluences. The 
equilibrium optical conductivity spectrum of the insulat- 
ing phase of VO2 is depicted in Fig. |6l[a) for comparison. 
Figs. IHb) and (c) show data for an excitation fluence 
$ below the transition threshold $c- As discussed in 
Sec. IIIIBl the threshold fluence depends on temperature 
Tl- In particular, excitation densities that completely 
switch the VO2 to the metallic phase at room tempera- 
ture are not sufficient to do so if the lattice is cooled down 
to 4 K. Performing the 2D measurement with a fluence 
of $ = 7.5 mJ/cm'^, at Tl = 4 K, results in the data 
displayed in Fig. [UJb). The observed dynamics are es- 
sentially identical to below-threshold excitation with the 
fluence $ = 3 mJ jcvc? at Tl = 250 K shown in Fig. Eljc) . 

Within our time resolution of 40 fs, ultrafast photo- 
doping induces a quasi-instantaneous onset of conductiv- 
ity due to directly injected mobile carriers. For pump 
fluences below the threshold, the electronic part of the 
optical conductivity Acri(aj,r) decays promptly within 
approximately 400 fs. In contrast, the pump-induced 
changes in the phonon resonances are more long-lived. 
While the onset of the phononic response varies for the 
three modes, we find intriguing common features. Pho- 
toexcitation induces increased polarizability on the low- 
frequency side of each phonon resonance while minimal 
change is seen on the blue wing: the resonance frequen- 
cies are red-shifted. For all modes, the change in fre- 
quency is superimposed on a remarkable coherent os- 
cillation of Acri(a;,T), along the pump-probe delay axis 
T. This phenomenon is most notable at a THz photon 
energy of 60 meV marked by a vertical broken line in 
Fig. [6l The corresponding cross sections of the 2D maps 
are shown in Fig [7l[a) . 

After subtraction of a slowly varying background, the 
coherent oscillation is well fitted by a cosine function with 
a frequency of 6.1 THz, as shown in Fig [Tljc). Thus, 
the true origin of the periodic conductivity modulation 
that appeared already in the spectrally integrated data 
of Fig. [5] is related to the periodic change in the phononic 
response. The analysis shows that the failure of the sim- 
ple decomposition in explaining the spectrally integrated 
conductivity dynamics for r < 400 fs [first oscillation cy- 
cle in Fig. EJa)] arises from the presence of the quickly 
decaying electronic contribution in the total response. 

As argued in Sec. IIII Bl the coherent oscillations at 
the single frequency of 6.1 THz observed by multi-THz 
probing in contrast to the all-optical experiments^ [Fig-E] 
indicate that the utilized multi-THz probe is more sen- 
sitive to the coherent dynamics of the excited electronic 



state of VO2. The 2D maps clearly demonstrate that 
the origin of the coherent dynamics traced by the multi- 
THz probe originate in the anharmonic coupling between 
the high-frequency oxygen and low-frequency vanadium 
lattice vibrations. Since the multi-THz response is dom- 
inated by the vanadium vibration inherent to the rutile 
phase, we assume that the anharmonicity of the VO2 lat- 
tice in the excited electronic state is strongly enhanced 
as compared to the insulating ground state. 

The spectral response of VO2 and its dynamics shown 
in Figs, inijb) and (c) are very similar. This fact once 
more convincingly demonstrates the interchangeability of 
temperature and fluence. We may thus conclude, that 
the observed interplay of electronic and phononic degrees 
of freedom that eventually leads to the emergence of the 
metallic phase, is not unique to the optically driven case 
but plays an important role in the thermal transition as 
well. 

Fig. El^d) shows the 2D map measured with an exci- 
tation fluence of $ = 7.5 mJ/cm^ exceeding a threshold 
fiuence $c = 5.3 niJ/cm'^ at a temperature Tl = 250 K. 
The spectral response of the optical phonon resonances is 
qualitatively similar to Figs. [Sfb) and (c) with the same 
dynamics characterized by the coherent oscillations [see 
Fig. El^a)]. In contrast, the dynamics of the electronic 
photoconductivity differs profoundly: After a resolution- 
limited onset of ACTi(a;,T) due to nearly instantaneous 
photodoping, the THz conductivity levels off at a high 
value, indicating the final transition into a metallic phase. 
This behavior is seen also in Fig. ^h) where the cross 
sections of the 2D maps at the energy of 100 meV are 
shown. 

Remarkably, at this energy there are no pronounced 
oscillations imprinted on the electronic conductivity, in 
contrast to the phononic response. Nevertheless, careful 
inspection reveals a periodic modulation of the multi- 
THz response around 90 meV. This energy range is well 
above the highest phonon resonance and related to the 
electronic part of the conductivity. Most clearly, this ef- 
fect is seen in the cross section of Fig. ^d) made at an 
energy of 92 meV. The oscillating component in this cross 
section shown in Fig. [Tjd) occurs exactly at twice the V- 
V stretching frequency of 6.1 THz. Moreover, the first 
three cycles of the oscillation shown in Fig. ^d) are well 
described by squaring the fitting function of the 6.1 THz 
oscillation in Fig. [Tjc) (and subtracting an offset). For 
later delay times r the modulation amplitude rapidly de- 
creases and a relation to the coherent wave packet mo- 
tion cannot be traced any more. The doubled frequency 
of the oscillation indicates that the dominant term in the 
electron-phonon coupling of the 6.1 THz phonon mode is 
quadratic with respect to the phonon normal coordinate. 

The red-shift of the oxygen modes is equal in all three 
investigated cases and hence depends neither on fluence 
nor temperature. The same holds for the period of the 
modulation, which always consistently matches the vana- 
dium 6.1 THz mode, as displayed in Fig. [Tja). Interest- 
ingly, the amplitude of the modulation scales with fluence 
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FIG. 8: (color online). Schematic energy surfaces as a func- 
tion of a dimer coordinate Q. The energy minima of the 
potential surfaces Qm and Qr correspond to spatial configu- 
rations of the V-V dimers in the monoclinic and rutile phases, 
respectively. Vertical dashed arrows denote Franck-Condon- 
type photoexcitation of spin singlets into a conductive state. 

(a) Structural relaxation and coherent vibrations about the 
new energy minimum below the pumping threshold $ < "I>c. 

(b) Intense excitation above the threshold $ > $c leads to a 
structurally assisted collapse of the mobility edge. 



but hardly depends on temperature at all. 



IV. DISCUSSION 

A qualitative understanding of the experimental ob- 
servations may be reached using a recent model that de- 
scribes the electronic structure in thermal equilibrium 
within a cluster dynamical mean- field theory.^" In this 
picture, the low-temperature insulating phase is regarded 
as a molecular solid of vanadium dimers embedded in a 
matrix of oxygen octahedra. To first order, the corre- 
lated electronic state of each dimer may be described by 
bonding and anti-bonding Heitler-London orbitals. The 
energy dependence of these states on the nuclear posi- 
tions (such as V-V separation) is depicted schematically 
in Fig. [51 The minimum of the lower-energy bonding 
state defines the atomic position in the monoclinic phase. 

Absorption of a near-infrared photon removes an elec- 
tron from the bonding orbital, destabilizing the dimer, 
while the lattice site is left in a vibrationally excited 
Franck-Condon state. Symmetry requires that the en- 
ergy minimum of anti-bonding orbitals should be located 
near the rutile configuration. Ultrafast photoexcitation 
thus launches a coherent structural deformation of ex- 
cited dimers followed by oscillations of Ag symmetry 
around the new potential minimum (Fig. [5]). The ap- 
preciable lattice deformations impose strong distortions 
on the surrounding oxygen octahedra and shift their res- 
onance frequencies. A similar effect has been reported in 
carbon nanotubes:^ Remarkably, the frequency shift is 
very pronounced for the oxygen vibrations polarized per- 
pendicular to the aj\/ axis (modes at 62 and 74 meV) and 
relatively weak for the phonon mode at 50 meV polarized 
along the aM axis. This fact attests to an anisotropic 
character of the anharmonic coupling between normal 



modes of the vanadium and oxygen sublattices. 

Since the oscillation period of the oxygen-related 
phonon modes is shorter by a factor of three as com- 
pared to the inverse frequency of the Ag vibration, the 
eigenfrequencies may follow the structural changes of the 
vanadium dimer adiabatically. Our experiment directly 
demonstrates the infiuence of a coherent lattice motion 
on other phonon resonances: An overall red-shift of the 
oxygen modes attests to a modified average structure 
of the vanadium dimers, while the coherent modulation 
at 6 THz refiects large-amplitude oscillations about the 
new potential minimum in the excited state. Although 
additional shifts of infrared-active phonon modes may 
be expected from screening via delocalized charge carri- 
ers, the distinctly different decay times of electronic and 
phononic conductivity prove that this effect is not domi- 
nant [Fig.[7] . The density of excited dimers scales linearly 
with the laser fluence and so does the amplitude of the 
observed coherent modulation. On the other hand, the 
lattice temperature has no measurable influence on the 
coherent modulation of the oxygen modes. 

The model also provides an instructive explanation of 
the transient electronic conductivity. Photodoping leaves 
the dimers initially in a higher-energy state with delo- 
calized electronic character (Fig. [5]) generating a quasi- 
instantaneous onset of electronic conductivity. Structural 
distortion drives the dimers into a new energy minimum. 
The rapid non-exponential decay of the electronic con- 
ductivity seen in Fig. [Tljb) indicates that the structural 
deformation shifts the energy of excited dimers below a 
mobility edge for extended electronic states [Fig. IDJa)]. 
This process is analogous to a self-trapping of excitons 
considered by Dexheimer et al»^ Although self-trapped 
states do not contribute to the electronic conductivity, 
their structural distortions manifest themselves by the 
anharmonic shift of oxygen-related phonon modes that 
persist for at least 1 ps [see Fig. [SJa)]. On the other 
hand, at high enough pump fluences electronic conductiv- 
ity does not decay rapidly. We suggest long-range coop- 
erative effects to account for this finding: Strong distor- 
tions that locally map onto the rutile lattice are expected 
to adiabatically renormalize the electronic bands^ and 
lower the mobility edge. Due to the enhanced critical 
fluctuations at elevated temperatures the system is more 
susceptible to this renormalization in accord with exper- 
imentally observed lowering of the critical fluence upon 
approaching T^. 

Instead of the persistent coherent modulation at the 
frequency of the Ag vibration, a weak modulation of the 
electronic conductivity at the doubled frequency is ob- 
served. The amplitude of this oscillation rapidly van- 
ishes after about three oscillation cycles [Fig.[SJd)]. This 
behavior suggests that the lowest-order term in the cou- 
pling between the Ag mode and the electronic conduc- 
tivity band of the photoexcited state is proportional to 
AQ^, where AQ = Q — Qr denotes the deviation of 
the dimer coordinate from the metastable rutile struc- 
ture (see Fig. |S]). Since, in a flrst approximation, this 



9 



coupling is rather weak, the electronic system may be 
considered as decoupled from the ionic motion in the 
non-equilibrium regime after the photoinduced insulator- 
metal transition. Consequently, the long-lived high value 
of electronic conductivity cannot originate from the con- 
tinuing motion of the V-V dimers. Instead, we propose 
to explain this feature by the collapse of the mobility 
edge^ below the Fermi energy [Fig. [Sfb)]. In this pic- 
ture, the critical pump fluence $c determines whether 
self-trapping or cooperative delocalization of electron- 
hole pairs prevails. 

V. CONCLUSION 

The ultrabroadband THz experiments presented here 
allow us to directly trace the the temporal evolution of 
the electronic conductivity and the optical phonon res- 
onances during the transient formation of the metallic 
phase in VO2. This study reveals fingerprints of a co- 
herent V-V intradimer wave packet motion at 6.1 THz 
which couples anharmonically to infrared-active phonon 
modes. The coherent response of the crystal lattice is ex- 
plained within a model description of local V-V dimers 
photoexcited into an antibonding state. The electronic 
conductivity shows a weak quadratic coupling to the co- 
herent wave packet motion. A threshold fluence of the 
excitation pulse for the photoinduced transition is found 



to decrease as the equilibrium lattice temperature Tl ap- 
proaches the transition temperature Tc- Below threshold, 
the electronic correlations remain undisturbed and the 
mid-infrared conductivity vanishes on a sub-picosecond 
time scale. Above a critical excitation density, a cooper- 
ative insulator-metal transition occurs on a timescale set 
by the first half of the oscillation cycle of the V-V coher- 
ent wavepacket motion in the excited state. In so far, the 
data indicate a strong structural component of the phase 
transition mechanism. However, right after the transition 
to the metallic phase the electronic system is rendered in- 
sensitive to the continuing coherent wave packet motion. 
Thus the stabilization of the metallic phase must origi- 
nate from a different mechanism, most likely a collapse 
of the mobility edge due to the cooperative modification 
of the electronic band structure. 
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